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One of the most exciting features of the hole centers CuOl~ in doped cuprates is an unusually 
complicated ground state which is the result of the electronic quasi-degeneracy. An additional hole, 
doped to the basic CuO^~ cluster with the big hole can occupy both the same hybrid Cu3d — 02p 
orbital state resulting in a Zhang- Rice singlet ^Aig and the purely oxygen Cu molecular orbital 
resulting in a singlet or triplet ^'"^Eu term with the close energies. We present detailed analysis of 
the (pseudo)-Jahn- Teller effect driven by the near-degeneracy within the ^Aig,^'^ _B„-manifold. 



■ I. INTRODUCTION. 

' Intuitive ideas concerning a specific role of Jahn- Teller ions (centers, polarons) have been used as a starting point 
J> of the pioneer investigations by K.A.MuUer and J.G. Bednorz resulting in 1986 in the outstanding discovery of the 
O high-Jc superconductivity [Q . Unfortunately, in the following years no breakthrough in understanding of this puzzling 
, phenomenon occurred. In many respects such situation is accounted for the underestimation of Jahn- Teller-effect and 
related phenomena that is typical for the conventional "metallic" approaches to the description of the electronic 
1*^^ , structure of the cuprates: namely these approaches underly the majority of the popular scenario's for the high-Tc 
superconductivity. Additional argumentation of the opponents of the Jahn- Teller approach is based on the widespread 
' opinion that the CuOa cluster with the doped hole forms a well-isolated spin and orbital singlet ^Aig (Zhang-Rice 
J> : singlet). 

CO ■ Well-isolated Zhang- Rice singlet is a natural starting point for many model approaches including the well-known t-J- 
, model. However, it does not provide an explanation of wide set of unconventional physical properties of the cuprates, 
' in particular, those associated with anomalous lattice and electron-lattice effects, and should be reconsidered and 
^ [ generalized. 

QQ ' A considerable number of experimental data were accumulated which more or less directly give evidence in favor 
of existence of the copper-oxygen centers with near-degeneracy effects and anomalously strong electron-vibrational 
. (or Jahn- Teller) correlations. More and more experimental data argue that a single band picture for the low energy 

■ excitations in the cuprates is inadequate. 

^ , An observation of the so-called mid-infrared (MIR) absorption bands in different CM04-cluster based oxides is one 
I of the most impressive manifestations of near-degeneracy effects |^ and, moreover, this phenomenon provides an 
' ^ important information about the energy spectrum and electronic structure within the ground state manifold. An 
' appearance of the MIR bands with specific transformation of the absorption spectra in cuprates upon doping leads 
^ ] to a conjecture that the latter is accompanied by the sharp correlational decrease in the energy of the charge transfer 
• • . transition big — > e„ that determines the fundamental absorption band for the parent oxides Thus, "parent" 
. ^ ' absorption band shifts from rather usual position around ^ 2 3 eV to the mid-infrared region forming the MIR 
, band. An appearance of the MIR bands upon hole doping could be readily explained by assuming that additional 
' hole doped to the basic CuO^ cluster with the big hole can occupy both the same hybrid Cu3d — 02p orbital state 
5^ , resulting in the Zhang-Rice singlet ^Aig and the purely oxygen e„ molecular orbital resulting in the singlet or triplet 
^'^Eu term with the close energies. Then the MIR-absorption is determined by the allowed charge transfer transition 
big : ^Aig — > bigCu : and represents the correlation analogue of the corresponding single-particle big —> Cu 

transition (see Figj^). 

Perhaps, a detection of isolated PJT center would provide the direct manifestation of the validity of the Jahn- 
Teller conception. In this connection the paper Q should be noted where the authors have performed NQR study 
of the isolated hole centers in La2Cuo.5Ljo.5O4. The results can be interpreted as convincing evidence of the quasi- 
degenerated singlet-triplet structure of the hole center. This conclusion is based on the following: 

1. The authors have detected spin singlet ground state {S = 0) and low lying spin triplet state {S — 1) with 
singlet-triplet separation Ast = 0.13 eF. 

2. They observed anomalously weak temperature dependence of the relaxation rate at low temperatures that gives 
evidence of occurrence of the spinless multiplct structure in the CuO^ cluster ground state. 

3. They've found considerable spin contribution to the low temperature relaxation indicating the simultaneous 
occurrence of the ground state multiplet structure, sufficiently low singlet-triplet separation and intrinsic singlet- 
triplet spin-orbital mixing. 

4. They observed the relaxation inequivalence of the various Cu sites, which is quite natural for the PJT centers 
in the conditions of the static PJT effect. 
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The Jahn- Teller hole centers like CuO^~ with singlet-triplet quasi-degeneracy within ground state have been 
observed by ESR spectroscopy in LaSrAli-xCux04 which is isostructural to La2-xSrxCuOi |^. Moreover, Yu. 
Yablokov et al. 1^ conjectured that doped hole in the copper-oxygen clusters occupies a purely oxygen a2g(7r), or 
^2ti(7r) like orbitals. An important indication to the 02p(7r) nature of doped holes and, hence, to the occurrence of 
near-degeneracy for configurations like h\g and higeu{i^) was obtained by Yoshinari ||^ and Martindale et al. |^ after 
the analysis of the ^^O Knight shift data and temperature-dependent anisotropy of the planar oxygen nuclear spin- 
lattice relaxation rate in Y Ba2Cu'iOQ+x, respectively. Analogous conclusion could be drawn out of the comparative 
analysis of the temperature behavior for the nuclear spin-lattice relaxation rate in La2-xSrxCuOi All this implies 
a complicated nature of the ground state manifold for the CuOi center with a significant mixing of the Zhang- Rice 
singlet and some other molecular term, which symmetry should be distinct from ^Aig. This conclusion conflicts with 
the widespread opinion regarding the well isolation of the Zhang-Rice singlet. 

An important argument in favor of vibronic nature for ground state of Cu04 clusters with the participation of e„- 
orbitals in the 123 system was obtained after analysis of experimental data on EFG (electric field gradient) tensor for 
different nuclei in 123 system |jl^: non-contradictory description of the data implies a considerable (10%!) difference 
in electron density for 0(2) and 0(3) oxygens. It is unlikely that this result could be obtained without PJT effect. 

These and many other results of resonance (ESR, NQR/NMR) experiments being precise local probes cast doubt 
on the validity of popular conceptions which are widely used as a starting point for the analysis of resonance and in 
broader sense for many other physical effects in cuprates. 

A considerable number of the experimental indications of Jahn- Teller (vibronic) effects is associated with observation 
of lattice instabilities, ferroelectric [|ll|, pyro- and piezoelectric anomalies, local static and dynamic distortions 
p3[ , various phonon anomalies and manifestation of substantial electron-phonon effects such as the line shift and 
Fano effect for the phonon modes associated with local PJT-active modes [^ |l^. These observations and many 
other lattice effects are signatures of unconventional strong electron-lattice interaction at work in the cuprates |l^] 
with highly nonlinear and nonadiabatic intrinsic dynamics. Up to now these phenomena are often considered as 
convincing evidence in favour of the " structural" scenario for the high-T^ superconductivity usually associated with 
polarons (pseudo- Jahn- Teller polarons) or bipolarons p^ , p^ . A number of the phonon anomalies could be associated 
with the vibronic pseudo-spin fluctuations or, in other words, with effects of the short-range cooperative Jahn- Teller 
ordering ||l^ . These anomalies are linked with the specific points both inside the BZ and on its boundary. 

In our opinion the indirect evidence for Jahn- Teller nature of the CuOi centers with active role namely of the copper- 
oxygen hybrid Qe,, mode was displayed by the maximum entropy method (MEM) in Y Ba2Cu-iOQ+x at x ^ 1 |l9| ]. The 
authors observed characteristic squarish deformation of the nuclear density for Cu atoms in the Cu02 plane due to the 
anomalously strong anharmonic low temperature {T = IbK) motion presumably of vibronic nature. Unconventional 
isotope effect and anomalous anisotropic pressure effect on Tc in doped cuprates could also be associated with vibronic 
effects. 

Thus, numerous experimental data show that Zhang- Rice model should be generalized with inclusion of near- 
degeneracy effects accompanied by PJT effect. 

Many authors have treated JT (or PJT) effect in doped cuprates in a rather general and various form both as a 
source of the local pairing and as a source of the unconventional physical properties ||2^-|2^. So, H. Kamimura 
p2[ proposed a mechanism of the HTSC due to the coherent bipolaron conduction induced by JT distortions. A 
dynamic Van Hove Jahn- Teller effect has been introduced by R.S. Markiewicz |^^. The essence of the vibronic model 
by M. Georgiev et al., L. Mihailov et al. is existence of off-centered apex oxygens as local polarizable states due 
to PJT effect. Their model is one of the numerous so-called anharmonic models of the HTSC based on the account of 
anharmonic motion of apical oxygen | p^ . However, it should be noted that recent studies by pulsed neutron scattering 
suggest that it is the in-plane Cu site rather than the apical oxygen site that may be split into two positions 
The absence of strong apical anomalies was stated earlier in Ref. . 

In the most cases JT (PJT) effect has been considered within standard e-problem for Cv?^ ion in octahedral 
environment (214 systems) or in square pyramidal environment pq] (123 systems). Instead of standard dx^-y^ , 
doublet some authors considered dx2-y2 , dxz, dyz manifold p7| |. However, all approaches originated from the 
assumption of near-degeneracy for some predominantly CwSd-states do not agree with experimentally observed large 
gap (^ 1.5 eV) separating ground dx2^y2 and any other CwSd-states. Besides, any JT model approach pretending 
to be universal should be originated from CuO^ center as the only common element of the crystalline and electronic 
structure of all the cuprates. 

A transformation of CuO^ clusters into PJT centers upon hole or electron doping to Cu02 planes is a principal 
element of the so called singlet-triplet PJT (ST-PJT) center model developed in P, p^ , p8| ~p2[ . In addition to cuprates 
like YBa2CuzOQ+x, La2~xSrxCuOi^ La2CuOi+s this model could be readily extended to a series of strongly corre- 
lated oxides like {K, Ba)Bi03, Lai^xSvxMnOs, La2NiOi+s including systems with the high- Tc superconductivity 
and colossal magnetoresistance. Their unconventional properties refiect a result of response of the system to non- 
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isovalent substitution that stabilizes phases providing the most effective screening of charge inhomogeneity. These 
phases in oxides can involve novel unconventional molecular cluster configurations like Jahn- Teller sp-center with 
anomalous high local polarizability and multi-mode behavior. 

The copper oxides based on CuOi clusters within this model are considered to be systems, which are unstable with 
regard to disproportionation reaction 

2CuOl~- ^ [CuOl-] + [CuOl-] (1) 

with formation of system of polar (hole - CuO\~ or electron - CuOl~) pseudo- Jahn- Teller (PJT) centers. These 
centers are distinguished by the so-called local S-boson or two electrons paired in the completely filled molecular 
orbital of the Cu04-cluster. In other words, the novel phase can be considered to be system of local spinless bosons 
moving in lattice of the hole PJT-centers [Cu04~] or the generalized quantum lattice bose-gas (or liquid) with 
boson concentration near Nb = 1/2. 

In a sense, this microscopic approach represents a particular generalization of Zhang- Rice model. 
Though without detailed analysis of the ^Aig,^'^ — aig — big — b2g — e„ vibronic problem, this model approach has 
been successfully applied for qualitative and semi-quantitative description of many physical properties of cuprates: 

331 , local structure distortions 



evelopment of the model and. 



MIR absorption bands B, isotope effect 1 32 1 , static and dynamic magnetic properties 
p^ , hyperfine coupling , phonon anomalies ]30| ], neutron scattering ]36| . Further d( 
first of all, possibilities of quantitative predictions implies a detailed analysis of the ^Aig,^'^ — aig — big — 629 — Gu 
vibronic problem. Moreover, this problem is of great independent importance as a non-trivial example of multi-mode 
PJT effect. 

Below, in our paper we consider in detail a vibronic structure of the isolated PJT center taking into account some 
effects associated with its singlet-triplet structure and spin-orbital coupling. First, in Section 2 a short consideration 
of the correlation driven near degeneracy effects for the Cu04-clusters will be done. Section 3 contains a detailed 
analysis of the ^Aig,^''^ E^ — aig — big — 629 ~ vibronic problem within CuOa cluster including an adiabatic 
potential, vibronic states and tunnelling effects for three different regimes. Some effects of spin-orbital coupling 
within Aig,^''^ E^) manifold are considered in Section 4. 



II. CORRELATIONS AND THE NEAR DEGENERACY EFFECTS FOR THE CUO4 CLUSTERS. 

At a glance the analysis of electronic structure and energy spectrum of the parent compounds such as 
La2-xMxCu04^^ Y Ba2C'u30e+x at a: = does not display any exotic peculiarities except quasi-two-dimensional 
antiferromagnetism determined by the strong exchange interaction for the big{dx2-y2) holes in the "basic" CuO\~ 
clusters. At the same time it is worth to pay attention to one important feature, namely to the exciton-band form of 
the fundamental absorption in the 1.5-^3.0 eV region strictly pronounced in the systems like R2C'uOi, Y Ba2Cu306+x, 
CuO (29). 

A peculiar character of this absorption connected with the allowed charge-transfer transition big — *■ e„ between the 
copper-oxygen 6ig-hybrid and the purely oxygen Cu-orbital provides evidence of the strongly correlated nature of the 
e„-electrons with formal occurrence of two types of the e„-states with and without strong correlation. This peculiarity 
is associated with the maximal hole density occurred for oxygen ions just in the e„-states of the Cu04-cluster and 
can be easily explained in the framework of the non-rigid anionic background model . This model introduces new 
correlation degree of freedom with two possible states of anionic background for the valent 02p-holes corresponding 
to two possible projections of the correlation pseudospin s = 1/2 and is described by simplified Hamiltonian 

Hcorr = Vi&x + Vsa^, (2) 

where Vi,3 are two electronic operators for the valent states. Simple approximation used in |37[ conjectures the linear 
n2p -dependence for ^1,3, where n2p is the 02p-hole number. According to optical data [^ , the correlation pseudospin 
sphtting can achieve the value ^ 0.5 eV. 

An increase in the 02p-hole concentration with hole doping CuOf^ — > CuO^~ results in a sharp increase of the Cu 
correlation splitting in the hole CuO|^-centers with transformation of the big — > fundamental band to the high- 
energy (big) {bigCu) subband and the low-energy big — > ^i^e^ subband to be well known as MIR (mid-infrared) 
absorption band 

Fig.|l| shows qualitatively the main results of the taking account of the considered "e„ correlations" for energy 
spectrum of basic CuO\~ cluster and hole CuO\~ cluster. It is important to mention an appearance in our model 
of two types of the orthogonal (!) molecular orbitals; for instance big (upper correlation sublevel) and b\g (lower 
correlation sublevel) states, e„ and e*. 
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Thus, we come to a conclusion about a near-degeneracy for two configurations b^g^ and 6*^6* with big and e* 
being the lowest correlation sublevels. This result does not drastically changes with taking account of an electrostatic 
interaction Vee between two holes. Moreover, just the Vee contribution was considered earlier |^ as a main reason 
for a near-degeneracy for the ^Aig and ^'^Eu terms formed by bl^ and ^igC* configurations. So, the both correlation 
efi'ects lead to a near-degeneracy in ground state of the hole center CuO\~ . 

Unusual properties of the (^^i^,^'^ manifold involving terms distinguished by the spin multiplicity, parity and 
orbital degeneracy provides unconventional behavior for the hole center CuO\~ with active interplay of various modes. 
As an extremely important result one should note that i?„-doublet has a nonquenched Izing like orbital moment that 
can be directed only along the C4-axis. 

A near degeneracy within Aig,^''^ E^) manifold can lead to conditions for the pseudo-Jahn- Teller effect with 
anomalously strong electron-lattice correlations with active local displacements modes of the (5e„ , Qbi and 
types. It should be noted here that the Qe^ modes are the only hybrid copper-oxygen modes, while the Qb^ and Qb^ 
modes are the pure oxygen ones. 



III. VIBRONIC COUPLING FOR THE Cf/O^ CENTERS. 



A. Adiabatic potential. 

Below we'll make use the notation |S'A/5r7) for basis wave functions from the {} Aig^^'^ E^) manifold (see Figj^). 
Here, S (= 0, 1), Ms are the total spin and its projection, r7 {— Aig, i?^, El^) labels the irreducible representation 
of symmetry group D^h for the CuOi center and its row, respectively, indicating the transformation properties of 
the orbital functions. Subsequently, we restrict ourselves with the linear vibronic coupling within the {"^ Aig,^'^ E^) 
manifold with the JT-active vibrational coordinates of the aig, big, b2g, e„ symmetry. 

Vibronic coupling for the isolated singlet or triplet ^'^E^ term has a well known for the E — bi — 52-problem p^ ] 
form diagonal in S and Ms 

Singlet terms ^E^ and ^ A\g interact due to linear vibronic coupling 

Joo^ig|v;,6|oo£;;) = ^ v,q,^ (4) 



determined by active vibrational coordinates Qej, Qe\- For the CuO^ cluster there are three normal coordinates with 
e„ symmetry, however, below we restrict ourselves with the choice of one active e„ vibration with an appropriate 
linear vibronic coupling constant Ve- 

The sum of elastic energy Vq , electronic Hamiltonian Vei , and vibronic Hamiltonian Vmb for the singlet S = spin 
manifold (^Aig,^Eu) with bare separation A^^ (see Fig.|l|) could be written as: 



2r,2 / -A + V,Q, V,Q, V,Qy 

UiQ)^y2^^-i+ \ VeQ. FaQa VpQ p , (5) 

i ^ \ VeQy VpQp -VMc ' 



where the indices both for the coupling coefficient and normal coodinates are defined as follows: aig ^ z, ^ x, 

el ^ big ^ a, b2g ^ /3, and A = Aae - V.qf. V, = y/^' - F/", ^ = V./^l 

An important information for the PJT center could be obtained with examination of the adiabatic potential (AP) 
surfaces e{Q), which are the roots of characteristic equation for U{Q). In our case this is reduced to a cubic equation 
with an extremely complicated expression for the roots. The coordinates of the minima Q^, energy and structure of 
electronic wave function at Q*^ and curvature of energy surface near can be obtained by the Opik-Prycc method 
p8[ , where an eigenvalue problem is treated only for extremal points of the AP. It should be noted that the Opik-Pryce 
method does not permit to find specific points of the AP without definite values of the derivative, however namely 
such a situation occurs for the upper sheets of the AP. The type of minima can be derived from curvature analysis 
for e{Q) near Q". 

The relations between the extremum coordinates and the wave function coefficients are 
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Q^ = -gfz^ Q2 = -2^)zx, Ql = -2qfhy, (6) 

where we have denoted the coefficients of decomposition of the electronic wave function as z, x and y for 
jOO^ig) , |00i?J) and |00i?|^), respectively. The eigenvalue problem for Vmb in the AP extrcmum points is given 
by: 

-A-2^5j,z2 -AE':,^zx -^Efjj^zy \ /z\ /z\ 

-4:E^jT zx -2E'}j, {x^ - y2) -AE% xy a; = A x , (7) 

-AE':jt zy -AE^jT xy 2EJj. {x^ ~ y^) ) \ V ) \V j 

where E'^jj, — \Viqf'' — /2ujf is the specific JT energy. 

The system complemented with a normalization condition x"^ + y'^ + z"^ — \ has 13 solutions, which are listed in 
the Table ||. Also the expressions for the quadratic form of the AP surface near the extremal points are given. These 
solutions could be divided into three groups. 

1. The first group NJT (non-JT) contains the only solution. Electronic part of wave function in the extremal point 
is pure jOOAig). The NJT-extremum is the minimum on lower sheet of AP, if 6 < (6 = —/S. + AEjrp — 2Ejrp) or on the 
upper sheet, if (A + 2Ejrp) < 0. In both cases the weak pseudo-Jahn- Teller effect takes place, when due to weakness 
of the vibronic coupling in comparison with the bare separation of electronic levels A, there are no low-symmetry 
distortions of the CuO^ cluster, and only the renormalization of the local e„ vibration frequency occurs: 



2. The second group JT^ (i — a, (3) contains four solutions, which are similar to results of the well-known E — bi — b2- 
problem. The wave function at the extremum points is a pure £"„ superposition. In further analysis the rhombic mode 
with larger JT energy will be called the "strong" (cr) mode, and that with smaller JT energy will be called the "weak" 
{a') one: Ejj, > Ejj,, {a, a' = a, /?); for the E — bi — 62-problem the Ejj, is the JT stabilization energy. Among four 
extremum points only the JT^ pair will correspond to minima. The minima are located on the lower sheet of the 
AP, if < (oo- = A -I- AEjj, - 2Ejj,), and on the middle one, if (A - 2Ejj.) > 0. The JT^/ pair represents saddle 
points. The wave functions at the minima are orthogonal each other. If E'^jj, = Ejj,, the equipotential continuum of 
minima, or trough, exists. The a mode frequency does not vary, and that of cr' mode is renormalized due to vibronic 
coupling: 

ojI^ujI, e=-^'(l-A.), X.-^. (9) 

It should be noted here that a type {Big or B2g) of the ground JT mode is of principal importance for the physics of 
the copper oxides. This is determined by the competition of vibronic parameters for the C'uSd — 02p and 02p — 02p 
bonds minimizing the Big and B2g modes, respectively. 

For the e„ vibrations the JT^ solutions correspond to the weak P JT effect: only renormalization of local e„ vibration 
occurs. If a mode is strong, the coordinates of minimum determine the rhombic distortion of the CuOa cluster along 
X- or y -direction (Fig.^). Accordingly, the softening of the or e*J vibration occurs; but as the both minima are 
equivalent, the frequencies of local modes remain twice degenerate. In a case of "strong" (3 mode with rectangular 
distortion of the CuOi cluster, the softening of that of e,j modes occurs, which co-directs to the cluster distortion. In 
all cases the expression for the renormalized local mode frequency is written as: 



3. The third group PJT^ (i = a, /3) includes eight solutions and corresponds to the most complicated case of the 
strong pseudo-Jahn- Teller effect. In this case the wave functions at the extremum points are the Aig — E^ hybrid 
states, and coefficients of superposition depend on the bare splitting A and the JT energies. The four from the eight 
PJTcr extrema (u is the "strong" rhombic mode) are minima, if Oo- > and 6 > 0. All the PJTcr minima are equivalent 
and allocated on the lower sheet of the AP. An arrangement of minima in space of the normal coordinates of the CuOi 
cluster is schematically shown in Fig.^ The wave functions at the minima are not orthogonal each other, that is a 
characteristic feature of strong pseudo-Jahn- Teller effect. The four PJTo-' extrema (a' is the "weak" rhombic mode) 
are saddle points. A specific case of degeneracy for the JT energies of rhombic modes should be examined separately. 
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The rhombic distortion of the CuOi cluster at the minimum points with the nonzero plane quadrupole moment 
is accompanied by the co-directed distortion with the electric dipole moment. The possible PJT^- distortions are 
shown in Fig.^j. 

Close to minimum M the equipotential surface of the quadratic form e{Q — is a five-dimensional ellipsoid 

with center located at Q^'''^^-' and principal axes to be turned with regard to basic ones. The a mode at minimum is 
mixed with co-directed e„ mode and aig mode giving rise to three local hybrid modes. The cr' mode is mixed with 
the second mode giving rise to two local hybrid modes. The mixing coefficients are proportional to the appropriate 
vibronic coupling constants. 

In a case, when Vz = 0, the frequencies of normal local hybrid modes are written as: 



Q.sQi): u;l=---^±J['-—^]+C^, (11) 




where 



A = Ujj{l~vl) , OjI, {1~ pI) , C = UJeUJa'VaPa, 



and 



(Qa,Q2) : 




(12) 



where 



E'^jrp a^b — b 

F/j, {a^ + bf' " aa + b' 



For the minima 1 and 3 Qi = Qx^ if cr = a (Qi, if cr = /3), Q2 = Qy, if cr = o; {Q2, if = /3); for minima 2 and 
4 Qi = Qy, if cr a {Q2, if cr = (3), Q2 = Qx, if cr = a {Qi, if cr = /3). Due to an equivalence of minima all local 
frequencies coincide. 

The type of minima on the lower AP sheet will be determined mainly by following quantities: 

a„ = A + 4F}r - 2E'}rp , b=-A + AEfjj. - 2^}^ (13) 

in the following way: 

a) NJT, if a,^ > and 6 < 0; 

b) JT^, if < and 6 > 0; 

c) PJT^, if > and 6 > 0; 

d) NJT and JT^, if a,, < and b < 0. 

The diagram of states of the lower AP sheet in space of parameters A, Ejj. and Ejj. at constant value of E^f^ is 
shown in Fig.^. A cross-section of the parameter space for constant value of E'jj, is shown in Fig.|^a. If A ^ +00, the 
lowest ^Aig level is well isolated and the lower AP sheet has a trivial NJT minimum. In contrast, at A ^ —00 the 
^Eu term becomes the lower one and usual E — bi — 62-problem with two JT^ minima on the lower AP sheet occurs. 
When coming together, the ^Aig and ^E^ terms are mixed by the e„ mode, and the e„ frequency is renormalized. 
This is accompanied by a formation of four PJTo- minima or three (NJT-t-JTo-) minima depending on magnitude of 
Ejj,. These two possibilities correspond to lines Ejj, — E^f^ and E'^j, — Ej'^ in Fig.^. 

With a motion along the line Ejj, = E'^jl^ from —00 up to +00 the curvature of the JTq- minima along e,i directions 
decreases up to zero on the line = Q. At this point the JTcr minima transform into saddle points with simultaneous 
appearence of the PJT^ minima (one of the JT^- minima splits along Qi, another one splits along Q2 direction). 
Moving from the line ao- = up to 6 = the magnitude of rhombic distortion decreases up to zero on the line b — Q. 
The magnitude of e„ coordinates of minima firstly increases, then decreases up to zero, reaching a maximum on the 
line A = —Ejj. + Ejj,. The Qz coordinate increases linearly (~ a„) from zero up to —qi^ on the line 6 = 0. Thus, 
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the four PJTct minima transform into one NJT minimum. Further, with A — > +00 the Aig and terms interact 
weaker and LUe — > i^e- 

The PJTct minima occur, if the e„ mode driven interaction of the Aig and terms is rather effective. When 
Ejj, > 'iEjj, — Ejj, (for example Ejrp = E'^jj,), the a mode driven interaction becomes more effective. In this case at 
6 = the NJT minimum appears on the fower AP sheet together with JT^ minima. All three minima have the same 
energy on the line A = —Ejrp + Ejrp. Further, at A +00 the JT^. minima become more flat, without varying their 
coordinates, and at Oo- = only the NJT minimum remains. 

In Figjsb the diagram of the upper AP sheets is shown. 



B. Tunnel splitting. 



The most complicated case of strong pseudo-Jahn- Teller effect can be treated in a framework of the tunnel Hamilto- 
nian p9[ ,when the localized vibrations at the AP minima are considered with taking account of the inter- well tunneling. 
It is supposed, that the minimum depth is larger than the typical phonon energy, hence the tunnel frequency is rather 
small. 



We consider the minimization problem for the total energy functional [5*] = ( 

H = fn + Vel + VQ + t>„6 , 



H 



5* ) with Hamiltonian 



(14) 



where Tq is kinetic energy of nuclei, Vei is electronic energy operator, Vq is elastic energy, Vyib is vibronic Hamiltonian. 
The is written as 



(15) 



M=l 



where Lpfj and x^m ^'^^ electronic (see Table and vibrational wave functions, respectively, centered at the minimum 



P JT^*^-* . The ground state vibrational wave function x^m ^ form 



Xm 



n 

fc=i 



,M 



exp 




M ST 




(16) 



where (w^'^)^ are the eigenvalues of a matrix of the quadratic form e {Q — Q'"'*^'), U^'^ is matrix for the unitary 
transformation to principal axes for e {Q — Q^'^'*^-'). If Vz — 0, the frequencies uj^ are equal to uj± Eq. (pT|), lo'^ Eq. 
( p^ ) and LOzj respectively. 

The variation of the energy functional E gives 



H c^E S c, 



(17) 



where H and S are the Hamiltonian Eq. (|lj) and overlap matrix, respectively. The solutions of the system Eq. ( |l7| ) 
give the tunnel states and the tunnel energy levels, respectively. 

Due to equivalence of the PJT^r minima the H and S matrices include only three types of the non-zero matrix 
elements: 

a) the diagonal matrix elements: 

H 



Xm 



Vm Xm I ~ \Vm Xm Wm Xm / - 



b) the non-diagonal matrix elements for the states with different both dipole and quadrupole moments: 



'Pm Xm' 



H 



M \ _ 



fM+l Xm+1 



9' 



Vm^ XmVm+1 Xm+i) = Sg; 



M \ _ 



c) the non-diagonal matrix elements for the states with different dipole moment but the same quadrupole one: 



"P m X m 



H 



^M+2 Xm+2) - ^d, \Vm Xm Wm+2 Xm+ 



Sd- 



The explicit expressions for these quantities at V"^ = are listed in Table ^ Thus, the system Eq. ( p^ ) is written 



as: 
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with eigenvectors 



H Hq Hd Hq \ 

Hq H Hq Hd 

Hd Hq H Hq 

Hq Hd Hq H ) 




E 



1 Sq Sd Sq 

Sq 1 Sq Sd 

Sd Sq 1 Sq 

Sq Sd Sq 1 




(18) 



Cl 
C4 



(1,1,1,1), Cl = -(1,-1,1,-1), 



1 

1 

V2 



(— sin 9, cos 0, sin 6, — cos ( 
(cos 0, sin 6, ~ cos 0, — sin ( 



(19) 



Vectors C3 and C4 are degenerated, hence there is a freedom in choice of 9. We assume = 0, then 



(20) 



Ca \Aig) XgU) + 



where functions | Eu ^'^, Eip'^^ coincide, respectively, with {|00i?,^) , |00i?{f)}, if = a or with {{\OOEi^ 
\00Ey))/V2, (- |00£;;) + \00Ey))/V2}, if o- = /?. Symmetric combinations of vibrational functions are: 



Xl + X2 + X3 + X4 



1 



Xs=2 



-Xi + X2 - X3 + X4 



Xei^) = ^ Ul 



X2 - X4 



X3 



(21) 



The symmetry E of the vibronic and vibrational functions coincide with that of a mode. With taking account of the 
normalization for tunnel states 



Nl,^{^Aj^Aj = l + 2Sq 

1-Sd 



Sd 



(22) 



Nl = = l-2Sq + Sd 

we come to following expressions for the tunnel energy levels: 



H + 2Hq+ Hd 



H-Hd 



l + 2Sq + Sd 



l-Sd 



H — 2Hq + Hd 

1 — 2Sq + Sd 



(23) 



These are shown in Fig.^ as a function of A with dimensionless coupling constant fee = 3. The frequencies of 
tunneling between the equivalent distorted configurations of the CuO^ cluster are determined by the splittings of the 
tunnel energy levels, which are much less than typical phonon energies. It should be noted, that for the pseudo-Jahn- 
Teller effect the symmetry of the ground vibronic and bare electronic states could be different unlike the conventional 
Jahn- Teller effect ("Ham's law") Fig.|| illustrates the situation, when A < (electronic ^Aig level is higher than 
^E.i^), but -Eyiig < -E-Bu- The reason is that in the pseudo-effect the vibronic interaction mixes different electronic 
levels, contrary to a case of the degenerated electronic states. It is worthy to note that in our case the vibronic E 
level is always higher in energy than the Aig and ones. 
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C. Degeneracy of JT energies for rhombic modes. 



When Ejrp — Ejj. (and with 0^ > 0, 6 > 0) the trough of minima appears on the lower AP sheet. It is a curve in 
four-dimensional space, which has following parametric form: 

= -<zi"'c,, = -2q(")c,d,cos|, Qj, = -2(7(°)c,d<,sin|, (24) 
Q„ = -2q^^^di COS./?, Qp = -2q^"^dl sirup, 

where = ^J^^^, da = \J'^^^^- The energy and wave function are 

£0 = -E-jT - \claa (25) 



2 



I*) = c„ 



\A,,)+da (cos||i?^)+sin||i?^)). 



respectively. The CuO^ cluster distortions as a function of ip are shown in Fig.^. 

For JTct case, when Ca- = 0, the situation is similar to the i? — e-problem with the doublet ground vibronic state. Its 
type are not changed in PJTg. case, if c^. ^ 1. However, with the increasing of from the one hand the continuum of 
E — e type splits into two parts in direction and from the other hand the Aig — E^ mixing reduces the symmetry of 
electronic wave function resulting in a singlet ground state. The potential energy in a small vicinity of the continuum 
minima is 

s = So + ^fi + ^ + ^-F.p, (26) 

where q^, r and p are count out from 

The radial vibration frequencies along the principal directions of the quadratic form Eq. (^6[) are and uj'_ Eq. 
(|l^), respectively. If ?'o + Po ^ ^/^'-^ £0 ^ then the energy and wave function of the ground vibronic singlet 

are 

2 (ro + r) (po + P) 

where xo is the ground state vibrational function. 

If the parameters approach to those typical for the NJT situation, the radius of trough is small, so the rotational 
term is not a small perturbation. The energy barrier in the higher symmetry point lowers, and the system turns into 
vibrational regime near the NJT minima with a singlet ground state, which is not described by the Eq. (p^). 

IV. VIBRONIC STATES IN PRESENCE OF THE SPIN-ORBIT COUPLING. 

Without taking account of the spin-orbit coupling the ^Eu term is isolated, and it has the six-fold degenerated 
ground vibronic state. The spin-orbit coupling mixes the Ms = states of the ^E^ and ^£'„ terms. As well, this 
coupling splits lower vibronic states of the ^E^ term, which have Ms = ±1. 

A. The Ms = states. 

1. Well-isolated ^'^Eu terms. 

If the ^Aig and ^'^Eu terms are well separated in energy {Aae ^ A_e, A), it is possible to consider the AP within a 
basis of the \00E^), \00E^), jlOi^J), \10Ey) states, and then take into account the vibronic coupling with the of ^Aig 
term as perturbation. The potential energy matrix U{Q) acquires a form: 
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/ -Ae + VMc 



\ 







-A 



E — VaQa 

-iX 






iX 

Ae + VaC 



-iX \ 


Ae - VcQa I 



(28) 



where A is submatrix element of the spin-orbit couphng. The eigenvalues and eigenvectors of tfiff) are written as 
follows 



£l=S- J(Ab+p)^ + A2, 



|1) = -zsinr?i |-p,0) + COS771 1/9, 1) , 
|2) = zsinr?2 |p,0) +cosr?2 1) , 



£2 = S - ^ (Ag - pY + A2 
es = T. + \J {Ae - pf + A^ , |3) = cosr/2 |p, 0) + isinr?2 1) , 
sa = ^ + \J {Ae + pf + X^ , |4) = cosr/i |-p,0) - isinr^i |p, 1) , 



(29) 



where 



(30) 



\p, S) = cos0\SOE^) + sme\SOEl) , 
\-p,S) =- sine \SOE^)+ cos e\SOEl), 

tan26' = tan2T?i = , 

VaQa -Ae - p 



tan2?72 



A 



-Ae + p 



The minima of the lower AP sheet are located on Qo- axis (with Ejj, > Ejj,, a, a' = a, l3) at points, which represent 
solutions of the equation: 



A- 



{V^\Q^\ + AeY 



(31) 



The Qr^y = is a point of discontinuity of the derivative. The non-trivial minima, which correspond to the low- 
symmetry cluster distortions, exist for the arbitrary large A, if A^ ^ 0, but with A ^ oo the minima depth and the 
distortions magnitude become negligibly small. The ground vibronic state is twice degenerated. In strong coupling 
scheme this is realized due to orthogonality of electronic states belonging to different minima of AP. In a case A <C Ejj, 
the expressions for the minimum points and their energy are written as follows: 



A2 



" ' ^ [ 2{Ae + 2E^j^)\ 



(32) 



£o = —Ae — Ejrp — 



2 {Ae + 2E^^) 

The mixing coefficient for the triplet spin states with the singlet ^Aig term wave function is proportional to 

A 



Ab + A. 



AE 



E' 



(33) 



JT 



2. Strong PJT-effect for the singlet spin states. 

In this case it is neccessary to consider the spin-orbit mixing of the timncl states with lower vibronic states of ^E^ 
term. The only non-zero matrix elements are those between the tunnel i^„-states and the lower vibronic states of ^E^ 
term with Ms = 0. The effective Hamiltonian matrix is 



/ -A'^ 






iX' 



-A% iX' 
-iX' 




-iX' \ 


^'e 
A'e J 



(34) 
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where 2A^ is an appropriate energy splitting, A' is a modified matrix element of the spin-orbit coupling: 

A'= AiV£l^/2d,(xi^^|x^'^)) , (35) 

where x'm X^o ^ the ground state vibrational functions, which are centered at the minimum M of PJT^ type 
and at the minimum of AP of the E ~bi ~ 62-problem with opposite to M sign of a mode, respectively. 



3. The singlet-triplet asymmetry of the vibronic coupling. 

Let the difference in the linear vibronic coupling constants for the singlet ^Eu and triplet ^Eu states corresponds to 
the following inequalities for the JT energies: E'}j,{^Eu) > Ejj,(^E^ and E'^r^i^E^) < E^j^{^Eu). Then the mmima 
of AP for the term arc located on the axis, and those of the '^E^^ therm are located on the axis. If 
the surfaces of AP for different terms intersect, than the taking account of the spin-orbit coupling could result in a 
complicated form of the AP with four minima. The states in minima located on the same axis are orthogonal to each 
other, and those located on different axis are not. Hence, even if A^; = the lower vibronic states are two doublets, 
which are separated by the tunnel splitting A. The frequency related to A corresponds to the combined pulsing 
motion of the electronic and nuclear density between h\g and 62^ distortions of the CuO^ cluster. 



B. The Ms = ±1 states. 

A joint operation of the vibronic and spin-orbital coupling for the Ms = ±1 states within the ^Eu term is described 
by the matrix 



VaQa^z + VnQpax + iMsXiay 



(36) 



where Ai is a submatrix element for the spin-orbital coupling within the orbital part of the ^E^ manifold, at are the 
Pauli matrices, and the energy is counted off the ^E^ manifold. 
The lower AP sheet has four extrema at points 



J0,±) 



= ±h^2Ey^{l-p1) , i = a,p, 



(37) 



where li = •^/fi./wj, Ejj, = Ejj,/hijJi, pi = Xi/2Ejrp, and g^"'^'' = at > 1. The parameter pi equals to the ratio 
of the Eu level splitting due to the spin-orbital coupling (2Ai) to that of due to the vibronic coupling (4£'j-p ). At 
bi| > 1 (« = a, /3) we come to the weak pseudo-Jahn- Teller effect with the only minimum at gj^'^^ = {i = a, f3) and 
effective local vibration frequencies 



(38) 



At < 1 ( Ejj, > Ejrp) the strong pseudo-Jahn- Teller effect occurs with the AP minima at q^''^\ and with the 

AP saddle points at g^?'^^ The effective frequencies for the local vibrations at the AP minima and corresponding 
energies are derived as follows: 



Col 



^l' (1 - K) , X„ 



E5t 



E' 



(39) 



s_ = -E^jT - hw^p. 



JT 



As is seen, imlike the usual E — bi — 62-problem, the a mode frequency is renormalized due to the spin-orbital coupling. 



The electronic wave functions ip^^^^ at the minima q^'^^ are reduced to the following form: 



r+ = 



(Ms) _ J_ 

V2 



(40) 
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In contrast to the above considered E — hi— 62-problem, the spin-orbital couphng now results in the non-orthogonality 
of the wave functions located in the minima: 



that leads to the tunnel splitting. 

Varying the energy functional with the basis functions 

where x± ai'e the vibrational functions centered at q^''^\ results in a matrix equation: 



H++ H+- 



E 



1 ^4 



where 



S 



(Ms) I (Ms) 



-i Ms Pa exp 



-2E-jr{l-pl)'' 



Pa 



2-Xa 



H+_ = {H++ - hujaB) S , B = E-jt{1- pI) (2 - pI) . 
The tunnel level energies and tunnel splitting are determined as follows: 



Eg = H++ - 



huJaB \S\ 



, Eu = H++ + 

B\S\ 



nujaB\s\ 



i-i^r 

Appropriate vibronic wave functions can be represented as: 



where 



In the limit |po- 1 <C 1 for the magnitude of tunnel splitting one obtains: 

Ai = 2Ai exp (-2E'}j,^ , 
that might be interpreted as a result of vibronic reduction for the purely spin-orbital splitting. 



(41) 



(42) 



(43) 



(44) 



(45) 



(46) 



(47) 



V. REDUCTION FACTORS. 



A concept of vibronic reduction is widely used in a theory of the Jahn- Teller effect. The reduction factor is equal 
to a ratio of the reduced matrix element of the electronic operator, which is calculated for vibronic ground states, to 
this one for the bare electronic states. In a case of the pseudo-Jahn- Teller effect the vibronic reduced matrix element 
is commonly a linear combination of the electronic ones. 

Within the strong coupling scheme the vibronic state is written as 

*r7 = ;^E^(rri) ^ ¥'ri7iXr.T2(ri7ir272|r7) , (48) 

Ti 7ir272 
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where r2 G F x Fi, Nr is the normalization factor; c(rri) is the amplitude of contribution to the vibronic function 
F, which is generated by the electronic Fi state; Xr272 is a symmetric linear combination of the vibrational states of 
equivalent minima; (ri7iF272|F7) is the Clebsh-Gordan coefficient. 
Consider the vibronic matrix element of electronic operator Vp- 



J2 <TT,) c{r'r[) 



(49) 



E(^rJ E (ri7i|f7r'i7y (ri7ir272|r7)(F'i7;r272|rV) 

r2 ii'('ii2 



Here we used the Wigner-Eckart-Kostcr 41 1 theorem for the electronic matrix element in a case of the simply 
reducible group: 



F7 



r7 



F'7' 



F 



F' 



F7|f7F'7' 



and the orthogonality relation for the vibrational functions: 

(xr272lxr^7^) = (Xr2) ^r2r;,^727^ 



(50) 



(51) 



where (xr2) normalization factor. Applying the Wigner-Eckart-Koster theorem for the vibronic matrix element 

we find for real representations: 



vl/p 



F F' 

Fi f; 



where 



T V 
Fi F'l 



(-1) 



r+r+r 



NrNr' 

E(-i)'Mx^2) 



ic(rri) c(r'F'i) x 



F F F' 
F'l F2 Fi 



(52) 



(53) 



Here [] is the 6F-symbol for point group G The factor (—1)'" is 1 for all irreducible represenation of D^h except 
(— 1)'^^3 = — 1 |43|. In the usual Jahn- Teller situation Eq. ( p3[ ) gives the well-known result for the vibronic reduction 
factors Eil. Note that 



F F' 
Fi F'l 



(_l)r+r'+ri+r; 



F' F 
F'l Fi 



(54) 



In our problem FijF'i — Aig,Eu\ F,F',F2,F'2 = Aig,I],i?„; only the operators of symmetry F = 
^igj ^2g, Big^ B2g, liavc the non-zero matrix elements. The non-trivial reduction factors are given in Table ^; the 
remaining ones could be derived with help of Eq. (|5^. The reduction factors as functions of A are shown in FigJ|a. 
One can see that E' operators are reduced very strongly. With an increase in the vibronic coupling the reduction 
factors tend to their limit values (Figj8|D), which can be obtained from Eq. ( p3| ) with iVr —> 1, {'X^^) 1- 

The components of external electric field parallel to the CuOi cluster plane {E _L C4) induce an electro-dipole 
transitions between tunnel Aig and E^ states. The reduced vibronic matrix element of the i?„ operator is written as: 



1. ixh) -> 



A 



19 



E,, 



(55) 



-((x^..) + (xk) 

1) and at Ca — d^, = l/\/2 the matrix element of electro-dipole 



Thus, at strong vibronic coupling (A'r 
transition is not renormalized. 

For the singlet Eu states an orbital contribution to the Zeeman energy is not zero only for the non-zero values of the 
^-component of external magnetic field. This situation could be examined with a help of an effective spin Hamiltonian 
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for the non-Kramers doublet 1^^. The vibronic coupling results in an essential renormalization of reduced matrix 
element of the ^-component of angular momentum: 



Eu 



Va, 



(56) 



In this case the vibronic reduction factor is proportional to overlap integral 5^' of the vibrational states centered at 
neighboring minima on the lower AP sheet. With the increasing in the vibronic coupling ^ 0, that leads to a 
complete quenching of the angular momentum. 



VI. SPIN HAMILTONIAN OF THE PJT CENTER. 



A detailed pattern of the energy spectrum in external magnetic field and magneto-resonance properties of the PJT 
centers are substantially depend on the bare A — E splitting A, vibronic parameters and relative magnitude of the 
vibronic and spin-orbital coupling. 

However, some common features are determined only by the symmetry considerations, in particular, by the C4 
axial symmetry and specific properties of the bare electronic basis functions. So, the orbital doublet ^'^Eu terms are 
characterized by the nonquenched highly anisotropic (Izing like) effective orbital moment 1 = 1/2 oriented only along 
the C4 axis. 

Without taking account of vibronic coupling the effective spin Hamiltonian for the ^'^E^ terms could be represented 
as follows: 

H = XilS,+f3H,gilJ + 2(3 ^ H,S,, (57) 

i—x^y^z 

where f3 is the Bohr magneton. Si are spin matrices (5* = for the spin singlet ^Eu term), gi is an effective orbital g 
-factor for the ^'^E^ term, which magnitude is determined by the structure of electronic function. Here, the first 
term describes the spin-orbital coupling, the second and the third ones correspond to the orbital and spin Zeeman 
coupling, respectively, with the purely isotropic spin g-tensor: gx = 9y = 9z = 2. Taking account of the spin-orbital 
^Eu-^Eu mixing leads to emergence of the spin anisotropy with additive contribution to the effective spin Hamiltonian: 

A2 



Van = DSl , D = ^st~ \J^l + « , (58) 

and to the effective axial anisotropy of the spin g-tensor: g^ = 2, g^ ~ gy = 2cos0, where 

A2 




Taking account of the vibronic coupling upon the conditions of weak pseudo-Jahn- Teller effect {pa- = Xi/2Ejrp > 1) 
does not vary a form of the effective spin Hamiltonian since the vibronic distortions are suppressed by the spin-orbital 
coupling. With the strong pseudo-Jahn- Teller effect ( = Xi/2Ejrp < 1) one should make use of the vibronic states 
(Eq), that results in a very complicated spin- vibronic effective Hamiltonian. A relatively simple situation occurs at 
p^-C 1, and small magnitude of the overlap for the vibrational functions, when, neglecting the spin-orbital ^E^-^E^ 
mixing, we come to an effective Hamiltonian: 

H = -PH.gi ix+lx-) U + /3 E ^^9^lSJ , (60) 

i—x.y,z 

where the overlap integral (x+|x-) is determined according to Eq. (Q) and for the axial g-tensor: g^ = 2 — Pa9i, 
9x = 9y = '^- It should be noted the principal difference in the (^-factor anisotropy with and without vibronic effects. 
The doublet of the states with different (±) projections of the effective orbital moment corresponds to the functions 
\ ^'1 ^''} and ^'g ^■'}, respectively, where, along with the above defined functions Eq. (^), we 

have introduced 

*f Xf+^\E^') XL°)), (61) 
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An upper label for the vibrational function underlines the difference in the location of the minima for the Ms = and 
Ms 7^ spin ^Eu states on the lower AP sheet. It should be noted also, that, in general, the overlap integral (x+|x_) 
depends on the z-component of magnetic field due to the orbital Zeeman coupling, which contributes to the energy of 
the lower sheet of the adiabatic potential. At PHz 3> Ai an expression for the corresponding energy acquires a form: 

e- = \ KQa + ^iQl) - ^{V^Qcf + iypQ^f + m^g^f. (62) 



VII. CONCLUSIONS 

We have presented a detailed analysis of the pseudo-Jahn- Teller effect within a bare electronic (^Aig,^''^ E^) manifold 
for the hole or electron CuOj^ centers in doped cuprates. Above we did not consider a number of problems generated 
by the occurrence of the PJT centers. Firstly, one should note strong inter-center coupling effects caused by the 
common oxygen ion and the related effects of the coupling with either tilting or buckling modes. The PJT centers 
are responsible for the numerous effects of the short-r ange o r long-range cooperative PJT ordering observed for the 
cuprates: some of them have been considered earlier [ pTpol . Note Ref. |^^, where some effects of the cooperative 
PJT ordering are considered within a model that could be readily modified for our scenario. Authors have performed 
a model calculation of the specific heat, the elastic moduli and the thermal expansion coefficient for the cuprates with 
cooperative PJT effect accompanied by the strong fluctuations of crystalline fields. 

An important problem is associated with an influence of the PJT centers to the local boson kinetics and super- 
conductivity. As a whole, this is an item for separate discussion though some effects of the vibronic reduction and 
isotope shift were briefly considered earlier js^ . 

In conclusion, we would like to pick up and shortly list again a number of experimental data which confirm namely 
an above developed specific scenario of the PJT centers in cuprates: 

1. Appearance of the MIR absorption bands for all considered cuprates. 

2. The numerous NQR-NMR reveal of the singlet-triplet near-degeneracy within ground state manifold for the hole 
centers. 

3. Observation of the anomalously strong anharmonic low temperature thermal motion of the copper atoms within 
the copper-oxygen hybrid Qe^ mode displayed by the maximum entropy method. 

4. Observation of the ferroelectric anomalies. 

5. Unusual copper isotope-shift effect in superconducting cuprates. 

6. Observation of the specific phonon anomalies. 

Among diverse peculiarities of the singlet-triplet PJT centers it should be especially emphasized a possible appear- 
ance of the so-called "tunnel paramagnetic centers" which spin states S = l,Ms = ±1 and S = IjMg = 0, 
respectively, are localized within different wells of the adiabatic potential. In other words, different spin states corre- 
spond to different local distortions of the CuO^ cluster. Spin dynamics and relaxation for the tunnel paramagnetic 
centers are cruiciably dependent on the magnitude and orientation of external magnetic field. These centers could be 
relatively readily transferred to the metastable state. An occurrence of the tunnel paramagnetic centers inside the 
small droplets of the PJT center phase was considered as an origin of the magnetization and magnetostriction 
anomalies in cuprate CuO. Perhaps, these are responsible for the unusual magnetic resonance signals observed in 
cuprate Eu2CuOi [Q. 

It should be noted that we did not undertake the task of reviewing all the available experimental data confirming 
the PJT nature of the CuO^ clusters in doped cuprates and comparing them with our model approach: it is a separate 
problem. At the same time, it should be noted that the PJT center model is entirely based on the vast amount of the 
experimental material. 

We consider the above results as an essential first step to the elaboration of the comprehensive theory of the PJT 
lattice in the doped cuprates. 
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TABLE I. 

Expressions of the decomposition coefficients for electronic wave function, extremum coordinates of AP surface and 
quadratic form of AP surface near the extremal points. The following notations are used: qf^^ =Vi/(jf,C(r = \f^^, 



da = ^ 9i — i'lx + 1v)/^i 92 = {—Ix + <ly)/V^- For the NJT extremum Eq = —A — Ejj,; for the JT^ extrema: 
Sq = -E'^jj,; for the PJTc^ extrema: Sq = -E'^j, - ^c^aa- 
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TABLE 11. 

The overlap and tunnel Hamiltonian matrix elements. The following notations are used; ctq 
70 = qi°hc„d^, 

qi = 7o cos ip — ao sin ^p, q2 = 70 sin tp + (Tq cos i/j, 
ill = cos^ + LV- sin^ ip, Q[ = cos^ "ip + sin^ V> 
O2 = W-i- sin^ ip + uj- cos^ if, = u}'_^ s\v? tp + u>'_ cos^ ip, 
CI3 = {—u}++ui-) sin (fi cos (fi, 03 = (—co'^ + ) sin ip cos f/', 
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TABLE III. The reduction factors: the f ; T, F'; Fi, F'l label the transformation properties of the electronic operator (with 
E = a at S' = /3, and E = /? at S' = a), the vibronic wave function and the electronic wave function, respectively. Additional 
notations: (x^,,) = 1 + + S^, (xl) = 1 - 25," + S^, (xL) = 1 - ■ 
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Figure captions 



FIG. 1. Correlation effects in the energy spectrum of the basic CuO^~ cluster and the hole CuO\~ center with numerical 
values (in eV) typical for oxides like CuO. On the right hand side we show a formation of the fundamental absorption spectra 
for the parent and hole doped oxides with peculiar MIR band in the latter case.The lower panel is an illustration of the bare (a) 
CuO'l^ center energy spectrum modification with taking into account the strong pseudo-.Talm- Teller-effect (b) and the tunnel 
splitting of the ground vibronic states (c). Distortions of the CuOa, cluster assotiated with the different AP minima are shown 
schematically in the insert. 

FIG. 2. Possible distortions of the CuOa cluster in the a) JTa minimum; b) JT^ minimum; c) PJTa minimum (the dipole 
moment is oriented along the CuO a cluster diagonal); b) PJT/3 minimum (the rotation angle value for the dipole moment with 
respect to that in PJTa minima is 4> = 7r/4). 

FIG. 3. The allocation of the PJT^ minima in a space of symmetric coordinates Qa,Qi,Q2 (c is the "strong" rhombic 
mode). If <T = a, then Qi = Qa=, Q2 = Qy; H cr = (3, then Qi = ^ = -Q.+Qy ^ 

FIG. 4. The diagram of the lower AP sheet states in the (A, -Ejt, -Ejt) space at certain fixed Ej^ = E"/^ value. The planes 
Oct = and 6 = divide the (A, Ejj^, -E}^) space into four parts corresponding to NJT-, JTo--, PJTo-- or NJT+JTo-type of the 
lower AP sheet minima. In the NJT+JT^ region the plane is shown, where the NJT- and JT^ minimum energies are equal. 

FIG. 5. a) The diagram of states of the lower AP sheet on the (A,_Ejy) plane at certain fixed _E}y and Ejj, value. The 
boundaries Ucr = and 6 = of the regions with different type of the lower AP sheet minima are shown. The lines Ejrp = -B'^^ 
and -EjT = ^jt correspond to the two different relative intensities of the vibronic coupling via a mode and Cu mode, b) The 
diagram of states of the upper AP sheets on the (A, Ejj.) plane at certain fixed E}j- and E}^ value. In the region to the left 
of the A = —2Ejrp line there is NJT minimum on the upper AP sheet and only the trivial non-analytical minimum on the 
middle AP sheet. In the region to the right of the A = 2Ejrp line there is JTa minima on the middle AP sheet and the trivial 
non-analytical minimum on the upper AP sheet. At -Bjt < 2i5}y the regions of parameters exist, where the PJTo- minima on 
the lower AP sheet coexist with the NJT minimum on the upper sheet or with the JT,, minima on the middle sheet. Between 
the A = —2E}j, and A = 2Ejrp lines only the trivial non-analytical minima on the upper sheets exist. 

FIG. 6. The A dependencies of a) the lower tunnel energy levels and b) the corresponding tunnel splittings at ke = 3, 
krr/ke = 0.5, k^i / k^r = 0.5, kz = where A;f = 2Ejrp/?kiJi. A significant increasing of the E level energy at the A interval border 
is connected with non-applicability of the chosen basis for a description of the exited vibronic states in the situation of shallow 
minima. 

FIG. 7. The possible CuOa cluster distortions in the equipotential continuum of minima as function of (p angle (see text). 
The n points correspond to the ion allocations ai ip — n-K 12. The a) and b) panels differ by the initial mutual orientation of 
the h\g and 629 distortions. The top and bottom fragments on each panel correspond to different initial mutual orientation of 

the {big,b2g) and distortions. The possible trajectories of the oxygen ions for the different values of 

shown. The radius of the copper ion trajectory is proportional to qi°^ . The equivalent CuO a cluster distortions with opposite 
phase of the copper ion {ncu — » ncu + 4) are also possible. 

FIG. 8. The reduction factor dependencies a) as a function of A at A;e = 3, k^/ke = 0.5, k^i/ka = 0.5, kz = Q 
[kf = 2Ejj./Tiuji); b) as a function of dimensionless coupling constant k at ke = k, ka/ke = 0.5, ka-i/ka = 0.5, kz = 0. 
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